Tutorial: GMOS Imaging Reduction - NGC 3310 - i-band

IMPORTANT: Edits are required for the tables below, and in STEP 1 and 2. After that, it should just work.

Datasets:

Science observation UT date(s) 2015 Dec 14
Data filename prefix(es) N20151214S

Target Frames

ScienceN20151214S239-242binning 2x2
Processed Biasfor 20151214N20151214S0346_bias.fits
Processed Twilight Flatfor 20151214N20151204S0446_flat_i.fits
Processed Fringe FrameN/AN/A

Fringe Frames

N/A

Twilight Flats

The twilight flat that matches the science target observations. Not enough on the 10th, so use some from the 4th too.

Processed: N20151204S0446_flat_i.fits
Twilight FramesN20151204S446-451binning 2x2
N20151210S359-361binning 2x2
Processed Biasfor 20151204N20151204S0501_bias.fits
Processed Biasfor 20151210N20151210S0430_bias.fits

Biases

The biases that match the science target observations and the twilight flats.

Processed: N20151214S0346_bias.fits
BiasesN20151214S346-350binning 2x2
Processed: N20151204S0501_bias.fits
BiasesN20151204S501-505binning 2x2
Processed: N20151210S0430_bias.fits
BiasesN20151210S430-434binning 2x2

STEP 0: Directory structure used by this Notebook

The directory paths used in this notebook are assumed to match this structure:

rootdir/
    calib/
        bias/
        flat/
        fringe/
    raw/
    redux/

STEP 1: Set up the Notebook and the PyRAF session

The data processing will happen in the Notebook. It is also convenient to have a PyRAF interactive session ready to use if needed, eg. to imexam. In this step, we will configure the Notebook and the PyRAF session.

This step needs to be repeated everytime the Notebook and/or the PyRAF session is restarted.

STEP 1.1: Set up the Notebook

STEP 1.1.1: Set up the variables.


In [5]:
import os
import os.path
import shutil

# WARNING: Make sure that the directory structure given in Step 0 
#          is already in place.

# EDIT HERE - Directories
rootdir = '/data/workspace/DRWorkshopAAO/NGC3310'

# EDIT HERE - File numbers and file names
biases = {
    'bias_20151214' : {
        'procbias' : 'N20151214S0346_bias.fits',
        'frames'   : [('N20151214S', '346-350'),
                     ]
    },
    'bias_20151204' : {
        'procbias'   : 'N20151204S0501_bias.fits',
        'frames'     : [('N20151204S', '501-505')],
    },
    'bias_20151210' : {
        'procbias'   : 'N20151210S0430_bias.fits',
        'frames'     : [('N20151210S', '430-434')],
    }

}

flats = {
    'flat_20151204' : {
        'procflat'   : 'N2015120402S0446_flat_i.fits',
        'frames'     : [('N20151204S', '446-451'),
                        ('N20151210S', '359-361')
                       ],
        'flatbias'   : [biases['bias_20151204']['procbias'],
                        biases['bias_20151210']['procbias'],
                       ] 
    }
}

#fringes = {
#    'fringe_20101111' : {
#        'procfringe'   : 'S20101111S0028_fringe_i.fits',
#        'frames'   : [('S20101111S', '28-30,40-43')],
#        'fringebias'   : biases['bias_20101111']['procbias'],
#        'fringeflat'   : flats['flat_20101102']['procflat'],
#    }
#}

science = {
    'procsci'   : 'ngc3310_i.fits',
    'sciroot'   : 'N20151214S',
    'sciframes' : '239-242',
    'scibias'   : biases['bias_20151214']['procbias'],
    'sciflat'   : flats['flat_20151204']['procflat'],
#    'scifringe' : fringes['fringe_20101111']['procfringe'],
}

logfile = 'NGC3310_i.log'

# Automatically set directories to match structure from Step 0
raw_path = os.path.join(rootdir, 'raw') + '/'
redux_path = os.path.join(rootdir, 'redux') + '/'
procbias_path = os.path.join(rootdir, 'calib', 'bias') + '/'
procflat_path = os.path.join(rootdir, 'calib', 'flat') + '/'
procfringe_path = os.path.join(rootdir, 'calib', 'fringe') + '/'
products_path = os.path.join(rootdir, 'products') + '/'

# Change to work directory
os.chdir(redux_path)

Step 1.1.2: Load and initialize the package

---> LAUNCH DS9 before proceeding further

Load the package required for the Notebook session and reset tasks to the default parameters.


In [2]:
from pyraf import iraf
iraf.gemini()
iraf.gmos()

iraf.unlearn(iraf.gemini, iraf.gmos, iraf.gemtools)

iraf.gmos.logfile = logfile

iraf.set(stdimage='imtgmos')


This is the initial release of the IRAF FITSUTIL package
to include support for FITS tile compression via 'fpack'.
Please send comments and questions to seaman@noao.edu.

 
     +------------------- Gemini IRAF Package -------------------+
     |             Version 1.13.1, December 7, 2015              |
     |             Requires IRAF v2.14.1 or greater              |
     |              Tested with Ureka IRAF v2.16                 |
     |             Gemini Observatory, Hilo, Hawaii              |
     |    Please use the help desk for submission of questions   |
     |  http://www.gemini.edu/sciops/helpdesk/helpdeskIndex.html |
     +-----------------------------------------------------------+
 
     Warning setting imtype=fits
     Warning setting use_new_imt=no
 
gemini/:
 f2/            gemtools/       gsaoi/          niri/
 flamingos/     gmos/           midir/          oscir/
 gemlocal/      gnirs/          nifs/           quirc/
gmos/:
 gbias          gfquick         gmosaic         gsappwave       gssdist
 gbpm           gfreduce        gmosexamples    gscalibrate     gsskysub
 gdark          gfresponse      gmosinfo        gscrmask        gsslitfunction
 gdisplay       gfscatsub       gmosinfoifu     gscrrej         gsstandard
 gfapsum        gfskysub        gmosinfoimag    gscut           gstransform
 gfcombine      gftransform     gmosinfospec    gsdrawslits     gswavelength
 gfcopy         giflat          gnscombine      gsextract       mostools/
 gfcube         gifringe        gnsdark         gsflat
 gfdisplay      giillum         gnsskysub       gsfquick
 gfextract      gireduce        gprepare        gsreduce
 gffindblocks   girmfringe      gqecorr         gsscatsub

STEP 1.2: Set up the PyRAF session

The commands shown here need to be run from the shell.

Open an xterm or a Terminal, then:

cd /data/workspace/DRWorkshopAAO/NGC3310/redux
pyraf

Then in the PyRAF session:

gemini
gmos
unlearn gemini
unlearn gmos
unlearn gemtools

Get your PyRAF configured. In the PyRAF session:

iraf.gmos.logfile = "NGC3310_i.log"
set stdimage=imtgmos

STEP 2: Create the file lists


In [6]:
os.chdir(redux_path)

# delete any pre-existing lists
iraf.delete('*.lis', verify='no')

# biases
for (root, frames) in biases['bias_20151214']['frames']:
    iraf.gemlist(root, frames, Stdout='tmp_'+root+'.lis')
iraf.concat('tmp_*.lis', 'bias_20151214.lis')
iraf.delete('tmp_*.lis', verify='no')

for (root, frames) in biases['bias_20151204']['frames']:
    iraf.gemlist(root, frames, Stdout='tmp_'+root+'.lis')
iraf.concat('tmp_*.lis', 'bias_20151204.lis')
iraf.delete('tmp_*.lis', verify='no')

for (root, frames) in biases['bias_20151210']['frames']:
    iraf.gemlist(root, frames, Stdout='tmp_'+root+'.lis')
iraf.concat('tmp_*.lis', 'bias_20151210.lis')
iraf.delete('tmp_*.lis', verify='no')


# flats
for (root, frames) in flats['flat_20151204']['frames']:
    iraf.gemlist(root, frames, Stdout='tmp_'+root+'.lis')
iraf.concat('tmp_*.lis', 'flat_20151204.lis')
iraf.delete('tmp_*.lis', verify='no')


# fringes
#for (root, frames) in fringes['fringe_20101111']['frames']:
#    iraf.gemlist(root, frames, Stdout='tmp_'+root+'.lis')
#iraf.concat('tmp_*.lis', 'fringe_20101111.lis')
#iraf.delete('tmp_*.lis', verify='no')


# science
iraf.gemlist(science['sciroot'], science['sciframes'], \
             Stdout='sci.lis')

# Get a full list
iraf.concat('*.lis', 'all.lis')

STEP 3: Visually inspect all the data


In [7]:
all_files = open('all.lis', 'r')
for line in all_files:
    image = line.strip() + '[1]'
    print image
    iraf.gdisplay(raw_path + image, 1, fl_paste='no')
    iraf.sleep(5)
    
all_files.close()


N20151204S0501[1]

GDISPLAY - Started: 2016-04-09T14:18:28

GDISPLAY: Temporary image name tmpout4365_69N20151204S0501.fits

GMULTIAMP - Started: 2016-04-09T14:18:28


GMULTIAMP - Finished: 2016-04-09T14:18:29
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:18:31

N20151204S0502[1]

GDISPLAY - Started: 2016-04-09T14:18:36

GDISPLAY: Temporary image name tmpout4365_90N20151204S0502.fits

GMULTIAMP - Started: 2016-04-09T14:18:36


GMULTIAMP - Finished: 2016-04-09T14:18:37
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:18:39

N20151204S0503[1]

GDISPLAY - Started: 2016-04-09T14:18:44

GDISPLAY: Temporary image name tmpout4365_111N20151204S0503.fits

GMULTIAMP - Started: 2016-04-09T14:18:45


GMULTIAMP - Finished: 2016-04-09T14:18:45
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:18:47

N20151204S0504[1]

GDISPLAY - Started: 2016-04-09T14:18:53

GDISPLAY: Temporary image name tmpout4365_132N20151204S0504.fits

GMULTIAMP - Started: 2016-04-09T14:18:53


GMULTIAMP - Finished: 2016-04-09T14:18:54
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:18:56

N20151204S0505[1]

GDISPLAY - Started: 2016-04-09T14:19:01

GDISPLAY: Temporary image name tmpout4365_153N20151204S0505.fits

GMULTIAMP - Started: 2016-04-09T14:19:01


GMULTIAMP - Finished: 2016-04-09T14:19:02
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:19:04

N20151210S0430[1]

GDISPLAY - Started: 2016-04-09T14:19:09

GDISPLAY: Temporary image name tmpout4365_174N20151210S0430.fits

GMULTIAMP - Started: 2016-04-09T14:19:09


GMULTIAMP - Finished: 2016-04-09T14:19:10
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:19:12

N20151210S0431[1]

GDISPLAY - Started: 2016-04-09T14:19:17

GDISPLAY: Temporary image name tmpout4365_195N20151210S0431.fits

GMULTIAMP - Started: 2016-04-09T14:19:17


GMULTIAMP - Finished: 2016-04-09T14:19:18
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:19:20

N20151210S0432[1]

GDISPLAY - Started: 2016-04-09T14:19:25

GDISPLAY: Temporary image name tmpout4365_216N20151210S0432.fits

GMULTIAMP - Started: 2016-04-09T14:19:25


GMULTIAMP - Finished: 2016-04-09T14:19:26
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:19:28

N20151210S0433[1]

GDISPLAY - Started: 2016-04-09T14:19:33

GDISPLAY: Temporary image name tmpout4365_237N20151210S0433.fits

GMULTIAMP - Started: 2016-04-09T14:19:33


GMULTIAMP - Finished: 2016-04-09T14:19:34
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:19:36

N20151210S0434[1]

GDISPLAY - Started: 2016-04-09T14:19:41

GDISPLAY: Temporary image name tmpout4365_258N20151210S0434.fits

GMULTIAMP - Started: 2016-04-09T14:19:41


GMULTIAMP - Finished: 2016-04-09T14:19:42
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:19:44

N20151214S0346[1]

GDISPLAY - Started: 2016-04-09T14:19:49

GDISPLAY: Temporary image name tmpout4365_279N20151214S0346.fits

GMULTIAMP - Started: 2016-04-09T14:19:50


GMULTIAMP - Finished: 2016-04-09T14:19:51
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:19:53

N20151214S0347[1]

GDISPLAY - Started: 2016-04-09T14:19:58

GDISPLAY: Temporary image name tmpout4365_300N20151214S0347.fits

GMULTIAMP - Started: 2016-04-09T14:19:58


GMULTIAMP - Finished: 2016-04-09T14:19:59
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:20:01

N20151214S0348[1]

GDISPLAY - Started: 2016-04-09T14:20:06

GDISPLAY: Temporary image name tmpout4365_321N20151214S0348.fits

GMULTIAMP - Started: 2016-04-09T14:20:06


GMULTIAMP - Finished: 2016-04-09T14:20:07
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:20:09

N20151214S0349[1]

GDISPLAY - Started: 2016-04-09T14:20:14

GDISPLAY: Temporary image name tmpout4365_342N20151214S0349.fits

GMULTIAMP - Started: 2016-04-09T14:20:14


GMULTIAMP - Finished: 2016-04-09T14:20:15
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:20:17

N20151214S0350[1]

GDISPLAY - Started: 2016-04-09T14:20:22

GDISPLAY: Temporary image name tmpout4365_363N20151214S0350.fits

GMULTIAMP - Started: 2016-04-09T14:20:22


GMULTIAMP - Finished: 2016-04-09T14:20:23
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:20:25

N20151204S0446[1]

GDISPLAY - Started: 2016-04-09T14:20:30

GDISPLAY: Temporary image name tmpout4365_384N20151204S0446.fits

GMULTIAMP - Started: 2016-04-09T14:20:30


GMULTIAMP - Finished: 2016-04-09T14:20:31
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:20:33

N20151204S0447[1]

GDISPLAY - Started: 2016-04-09T14:20:38

GDISPLAY: Temporary image name tmpout4365_405N20151204S0447.fits

GMULTIAMP - Started: 2016-04-09T14:20:38


GMULTIAMP - Finished: 2016-04-09T14:20:39
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:20:42

N20151204S0448[1]

GDISPLAY - Started: 2016-04-09T14:20:47

GDISPLAY: Temporary image name tmpout4365_426N20151204S0448.fits

GMULTIAMP - Started: 2016-04-09T14:20:47


GMULTIAMP - Finished: 2016-04-09T14:20:48
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:20:50

N20151204S0449[1]

GDISPLAY - Started: 2016-04-09T14:20:55

GDISPLAY: Temporary image name tmpout4365_447N20151204S0449.fits

GMULTIAMP - Started: 2016-04-09T14:20:55


GMULTIAMP - Finished: 2016-04-09T14:20:56
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:20:58

N20151204S0450[1]

GDISPLAY - Started: 2016-04-09T14:21:03

GDISPLAY: Temporary image name tmpout4365_468N20151204S0450.fits

GMULTIAMP - Started: 2016-04-09T14:21:03


GMULTIAMP - Finished: 2016-04-09T14:21:04
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:21:06

N20151204S0451[1]

GDISPLAY - Started: 2016-04-09T14:21:11

GDISPLAY: Temporary image name tmpout4365_489N20151204S0451.fits

GMULTIAMP - Started: 2016-04-09T14:21:11


GMULTIAMP - Finished: 2016-04-09T14:21:12
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:21:14

N20151210S0359[1]

GDISPLAY - Started: 2016-04-09T14:21:19

GDISPLAY: Temporary image name tmpout4365_510N20151210S0359.fits

GMULTIAMP - Started: 2016-04-09T14:21:19


GMULTIAMP - Finished: 2016-04-09T14:21:20
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:21:22

N20151210S0360[1]

GDISPLAY - Started: 2016-04-09T14:21:27

GDISPLAY: Temporary image name tmpout4365_531N20151210S0360.fits

GMULTIAMP - Started: 2016-04-09T14:21:27


GMULTIAMP - Finished: 2016-04-09T14:21:28
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:21:30

N20151210S0361[1]

GDISPLAY - Started: 2016-04-09T14:21:36

GDISPLAY: Temporary image name tmpout4365_552N20151210S0361.fits

GMULTIAMP - Started: 2016-04-09T14:21:36


GMULTIAMP - Finished: 2016-04-09T14:21:37
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:21:39

N20151214S0239[1]

GDISPLAY - Started: 2016-04-09T14:21:44

GDISPLAY: Temporary image name tmpout4365_573N20151214S0239.fits

GMULTIAMP - Started: 2016-04-09T14:21:44


GMULTIAMP - Finished: 2016-04-09T14:21:45
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:21:47

N20151214S0240[1]

GDISPLAY - Started: 2016-04-09T14:21:52

GDISPLAY: Temporary image name tmpout4365_594N20151214S0240.fits

GMULTIAMP - Started: 2016-04-09T14:21:52


GMULTIAMP - Finished: 2016-04-09T14:21:53
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:21:55

N20151214S0241[1]

GDISPLAY - Started: 2016-04-09T14:22:00

GDISPLAY: Temporary image name tmpout4365_615N20151214S0241.fits

GMULTIAMP - Started: 2016-04-09T14:22:00


GMULTIAMP - Finished: 2016-04-09T14:22:01
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:22:03

N20151214S0242[1]

GDISPLAY - Started: 2016-04-09T14:22:08

GDISPLAY: Temporary image name tmpout4365_636N20151214S0242.fits

GMULTIAMP - Started: 2016-04-09T14:22:08


GMULTIAMP - Finished: 2016-04-09T14:22:09
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY - Finished: 2016-04-09T14:22:11

STEP 4: Process the biases

This will produce the two biases defined above. Those two biases will be used later to process the science frames, the twilight flats, and the fringe frames.


In [8]:
os.chdir(redux_path)

for biaskey in biases.keys():
    iraf.imdelete(biases[biaskey]['procbias'], verify='no')
    iraf.imdelete('g//@'+biaskey+'.lis', verify='no')

    print 'Creating bias ', biases[biaskey]['procbias']
    iraf.gbias('@'+biaskey+'.lis', outbias=biases[biaskey]['procbias'], \
               rawpath=raw_path, fl_over='yes', fl_trim='yes', \
               fl_vardq='yes')
    shutil.copy(biases[biaskey]['procbias'], procbias_path)

# The processed biases are:
os.listdir(procbias_path)


Creating bias  N20151210S0430_bias.fits
----------------------------------------------------------------------------
GBIAS -- Sat Apr  9 14:22:22 HST 2016

Input images:
  N20151210S0430
  N20151210S0431
  N20151210S0432
  N20151210S0433
  N20151210S0434

Output bias image: N20151210S0430_bias.fits
 

----------------------------------------------------------------------------
GIREDUCE -- Sat Apr  9 14:22:23 HST 2016

----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:22:23 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0430
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0430.fits   Output  gN20151210S0430.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:22:24 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0431
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0431.fits   Output  gN20151210S0431.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:22:25 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0432
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0432.fits   Output  gN20151210S0432.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:22:27 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0433
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0433.fits   Output  gN20151210S0433.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:22:28 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0434
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0434.fits   Output  gN20151210S0434.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------

Input files:
gN20151210S0430
gN20151210S0431
gN20151210S0432
gN20151210S0433
gN20151210S0434

Output files:
tmpfile4365_662N20151210S0430
tmpfile4365_662N20151210S0431
tmpfile4365_662N20151210S0432
tmpfile4365_662N20151210S0433
tmpfile4365_662N20151210S0434

GIREDUCE: Image tmpfile4365_662N20151210S0430 overscan subtracted
GIREDUCE: Image tmpfile4365_662N20151210S0431 overscan subtracted
GIREDUCE: Image tmpfile4365_662N20151210S0432 overscan subtracted
GIREDUCE: Image tmpfile4365_662N20151210S0433 overscan subtracted
GIREDUCE: Image tmpfile4365_662N20151210S0434 overscan subtracted
GIREDUCE: Image tmpfile4365_662N20151210S0430 trimmed
GIREDUCE: Image tmpfile4365_662N20151210S0431 trimmed
GIREDUCE: Image tmpfile4365_662N20151210S0432 trimmed
GIREDUCE: Image tmpfile4365_662N20151210S0433 trimmed
GIREDUCE: Image tmpfile4365_662N20151210S0434 trimmed


GIREDUCE - Cleaning up -- Sat Apr  9 14:23:17 HST 2016

GIREDUCE -- Sat Apr  9 14:23:17 HST 2016
GIREDUCE exit status:  good.
----------------------------------------------------------------------------
--------------------------------------------------------------------------------
GEMCOMBINE -- Sat Apr  9 14:23:17 HST 2016
 
GEMCOMBINE: input = @tmpcombin4365_659
GEMCOMBINE: output = N20151210S0430_bias.fits
GEMCOMBINE: title = Bias
GEMCOMBINE: bpmfile = 
GEMCOMBINE: sci_ext = SCI
GEMCOMBINE: var_ext = VAR
GEMCOMBINE: dq_ext = DQ
 
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:23: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:23
                Images     Exp
  tmpfile4365_662N20151210S0430.fits[SCI,1]    0.0
  tmpfile4365_662N20151210S0431.fits[SCI,1]    0.0
  tmpfile4365_662N20151210S0432.fits[SCI,1]    0.0
  tmpfile4365_662N20151210S0433.fits[SCI,1]    0.0
  tmpfile4365_662N20151210S0434.fits[SCI,1]    0.0

  Output image = tmpcombout4365_727, ncombine = 5
  Number rejected mask = tmpcombdq4365_728.pl
  Sigma image = tmpcombsig4365_729
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:23: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:23
                Images     Exp
  tmpfile4365_662N20151210S0430.fits[SCI,2]    0.0
  tmpfile4365_662N20151210S0431.fits[SCI,2]    0.0
  tmpfile4365_662N20151210S0432.fits[SCI,2]    0.0
  tmpfile4365_662N20151210S0433.fits[SCI,2]    0.0
  tmpfile4365_662N20151210S0434.fits[SCI,2]    0.0

  Output image = tmpcombout4365_734, ncombine = 5
  Number rejected mask = tmpcombdq4365_735.pl
  Sigma image = tmpcombsig4365_736
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:23: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:23
                Images     Exp
  tmpfile4365_662N20151210S0430.fits[SCI,3]    0.0
  tmpfile4365_662N20151210S0431.fits[SCI,3]    0.0
  tmpfile4365_662N20151210S0432.fits[SCI,3]    0.0
  tmpfile4365_662N20151210S0433.fits[SCI,3]    0.0
  tmpfile4365_662N20151210S0434.fits[SCI,3]    0.0

  Output image = tmpcombout4365_738, ncombine = 5
  Number rejected mask = tmpcombdq4365_739.pl
  Sigma image = tmpcombsig4365_740
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:23: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:23
                Images     Exp
  tmpfile4365_662N20151210S0430.fits[SCI,4]    0.0
  tmpfile4365_662N20151210S0431.fits[SCI,4]    0.0
  tmpfile4365_662N20151210S0432.fits[SCI,4]    0.0
  tmpfile4365_662N20151210S0433.fits[SCI,4]    0.0
  tmpfile4365_662N20151210S0434.fits[SCI,4]    0.0

  Output image = tmpcombout4365_742, ncombine = 5
  Number rejected mask = tmpcombdq4365_743.pl
  Sigma image = tmpcombsig4365_744
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:23: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:23
                Images     Exp
  tmpfile4365_662N20151210S0430.fits[SCI,5]    0.0
  tmpfile4365_662N20151210S0431.fits[SCI,5]    0.0
  tmpfile4365_662N20151210S0432.fits[SCI,5]    0.0
  tmpfile4365_662N20151210S0433.fits[SCI,5]    0.0
  tmpfile4365_662N20151210S0434.fits[SCI,5]    0.0

  Output image = tmpcombout4365_746, ncombine = 5
  Number rejected mask = tmpcombdq4365_747.pl
  Sigma image = tmpcombsig4365_748
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:23: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:23
                Images     Exp
  tmpfile4365_662N20151210S0430.fits[SCI,6]    0.0
  tmpfile4365_662N20151210S0431.fits[SCI,6]    0.0
  tmpfile4365_662N20151210S0432.fits[SCI,6]    0.0
  tmpfile4365_662N20151210S0433.fits[SCI,6]    0.0
  tmpfile4365_662N20151210S0434.fits[SCI,6]    0.0

  Output image = tmpcombout4365_750, ncombine = 5
  Number rejected mask = tmpcombdq4365_751.pl
  Sigma image = tmpcombsig4365_752
 
GEMCOMBINE -- Finished: Sat Apr  9 14:23:29 HST 2016
 
GEMCOMBINE -- Exit status: GOOD
--------------------------------------------------------------------------------

GBIAS: Exit status: GOOD.
----------------------------------------------------------------------------
Creating bias  N20151214S0346_bias.fits
----------------------------------------------------------------------------
GBIAS -- Sat Apr  9 14:23:30 HST 2016

Input images:
  N20151214S0346
  N20151214S0347
  N20151214S0348
  N20151214S0349
  N20151214S0350

Output bias image: N20151214S0346_bias.fits
 

----------------------------------------------------------------------------
GIREDUCE -- Sat Apr  9 14:23:30 HST 2016

----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:23:30 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0346
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0346.fits   Output  gN20151214S0346.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:23:32 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0347
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0347.fits   Output  gN20151214S0347.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:23:33 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0348
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0348.fits   Output  gN20151214S0348.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:23:34 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0349
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0349.fits   Output  gN20151214S0349.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:23:36 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0350
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0350.fits   Output  gN20151214S0350.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------

Input files:
gN20151214S0346
gN20151214S0347
gN20151214S0348
gN20151214S0349
gN20151214S0350

Output files:
tmpfile4365_760N20151214S0346
tmpfile4365_760N20151214S0347
tmpfile4365_760N20151214S0348
tmpfile4365_760N20151214S0349
tmpfile4365_760N20151214S0350

GIREDUCE: Image tmpfile4365_760N20151214S0346 overscan subtracted
GIREDUCE: Image tmpfile4365_760N20151214S0347 overscan subtracted
GIREDUCE: Image tmpfile4365_760N20151214S0348 overscan subtracted
GIREDUCE: Image tmpfile4365_760N20151214S0349 overscan subtracted
GIREDUCE: Image tmpfile4365_760N20151214S0350 overscan subtracted
GIREDUCE: Image tmpfile4365_760N20151214S0346 trimmed
GIREDUCE: Image tmpfile4365_760N20151214S0347 trimmed
GIREDUCE: Image tmpfile4365_760N20151214S0348 trimmed
GIREDUCE: Image tmpfile4365_760N20151214S0349 trimmed
GIREDUCE: Image tmpfile4365_760N20151214S0350 trimmed


GIREDUCE - Cleaning up -- Sat Apr  9 14:24:26 HST 2016

GIREDUCE -- Sat Apr  9 14:24:26 HST 2016
GIREDUCE exit status:  good.
----------------------------------------------------------------------------
--------------------------------------------------------------------------------
GEMCOMBINE -- Sat Apr  9 14:24:26 HST 2016
 
GEMCOMBINE: input = @tmpcombin4365_757
GEMCOMBINE: output = N20151214S0346_bias.fits
GEMCOMBINE: title = Bias
GEMCOMBINE: bpmfile = 
GEMCOMBINE: sci_ext = SCI
GEMCOMBINE: var_ext = VAR
GEMCOMBINE: dq_ext = DQ
 
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:24: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:24
                Images     Exp
  tmpfile4365_760N20151214S0346.fits[SCI,1]    0.0
  tmpfile4365_760N20151214S0347.fits[SCI,1]    0.0
  tmpfile4365_760N20151214S0348.fits[SCI,1]    0.0
  tmpfile4365_760N20151214S0349.fits[SCI,1]    0.0
  tmpfile4365_760N20151214S0350.fits[SCI,1]    0.0

  Output image = tmpcombout4365_825, ncombine = 5
  Number rejected mask = tmpcombdq4365_826.pl
  Sigma image = tmpcombsig4365_827
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:24: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:24
                Images     Exp
  tmpfile4365_760N20151214S0346.fits[SCI,2]    0.0
  tmpfile4365_760N20151214S0347.fits[SCI,2]    0.0
  tmpfile4365_760N20151214S0348.fits[SCI,2]    0.0
  tmpfile4365_760N20151214S0349.fits[SCI,2]    0.0
  tmpfile4365_760N20151214S0350.fits[SCI,2]    0.0

  Output image = tmpcombout4365_832, ncombine = 5
  Number rejected mask = tmpcombdq4365_833.pl
  Sigma image = tmpcombsig4365_834
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:24: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:24
                Images     Exp
  tmpfile4365_760N20151214S0346.fits[SCI,3]    0.0
  tmpfile4365_760N20151214S0347.fits[SCI,3]    0.0
  tmpfile4365_760N20151214S0348.fits[SCI,3]    0.0
  tmpfile4365_760N20151214S0349.fits[SCI,3]    0.0
  tmpfile4365_760N20151214S0350.fits[SCI,3]    0.0

  Output image = tmpcombout4365_836, ncombine = 5
  Number rejected mask = tmpcombdq4365_837.pl
  Sigma image = tmpcombsig4365_838
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:24: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:24
                Images     Exp
  tmpfile4365_760N20151214S0346.fits[SCI,4]    0.0
  tmpfile4365_760N20151214S0347.fits[SCI,4]    0.0
  tmpfile4365_760N20151214S0348.fits[SCI,4]    0.0
  tmpfile4365_760N20151214S0349.fits[SCI,4]    0.0
  tmpfile4365_760N20151214S0350.fits[SCI,4]    0.0

  Output image = tmpcombout4365_840, ncombine = 5
  Number rejected mask = tmpcombdq4365_841.pl
  Sigma image = tmpcombsig4365_842
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:24: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:24
                Images     Exp
  tmpfile4365_760N20151214S0346.fits[SCI,5]    0.0
  tmpfile4365_760N20151214S0347.fits[SCI,5]    0.0
  tmpfile4365_760N20151214S0348.fits[SCI,5]    0.0
  tmpfile4365_760N20151214S0349.fits[SCI,5]    0.0
  tmpfile4365_760N20151214S0350.fits[SCI,5]    0.0

  Output image = tmpcombout4365_844, ncombine = 5
  Number rejected mask = tmpcombdq4365_845.pl
  Sigma image = tmpcombsig4365_846
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:24: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:24
                Images     Exp
  tmpfile4365_760N20151214S0346.fits[SCI,6]    0.0
  tmpfile4365_760N20151214S0347.fits[SCI,6]    0.0
  tmpfile4365_760N20151214S0348.fits[SCI,6]    0.0
  tmpfile4365_760N20151214S0349.fits[SCI,6]    0.0
  tmpfile4365_760N20151214S0350.fits[SCI,6]    0.0

  Output image = tmpcombout4365_848, ncombine = 5
  Number rejected mask = tmpcombdq4365_849.pl
  Sigma image = tmpcombsig4365_850
 
GEMCOMBINE -- Finished: Sat Apr  9 14:24:38 HST 2016
 
GEMCOMBINE -- Exit status: GOOD
--------------------------------------------------------------------------------

GBIAS: Exit status: GOOD.
----------------------------------------------------------------------------
Creating bias  N20151204S0501_bias.fits
----------------------------------------------------------------------------
GBIAS -- Sat Apr  9 14:24:39 HST 2016

Input images:
  N20151204S0501
  N20151204S0502
  N20151204S0503
  N20151204S0504
  N20151204S0505

Output bias image: N20151204S0501_bias.fits
 

----------------------------------------------------------------------------
GIREDUCE -- Sat Apr  9 14:24:39 HST 2016

----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:24:39 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0501
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0501.fits   Output  gN20151204S0501.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:24:41 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0502
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0502.fits   Output  gN20151204S0502.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:24:42 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0503
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0503.fits   Output  gN20151204S0503.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:24:43 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0504
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0504.fits   Output  gN20151204S0504.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:24:45 HST 2016

Input list    = /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0505
Output list   = 
Output prefix = g
Raw path      = 
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0505.fits   Output  gN20151204S0505.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------

Input files:
gN20151204S0501
gN20151204S0502
gN20151204S0503
gN20151204S0504
gN20151204S0505

Output files:
tmpfile4365_858N20151204S0501
tmpfile4365_858N20151204S0502
tmpfile4365_858N20151204S0503
tmpfile4365_858N20151204S0504
tmpfile4365_858N20151204S0505

GIREDUCE: Image tmpfile4365_858N20151204S0501 overscan subtracted
GIREDUCE: Image tmpfile4365_858N20151204S0502 overscan subtracted
GIREDUCE: Image tmpfile4365_858N20151204S0503 overscan subtracted
GIREDUCE: Image tmpfile4365_858N20151204S0504 overscan subtracted
GIREDUCE: Image tmpfile4365_858N20151204S0505 overscan subtracted
GIREDUCE: Image tmpfile4365_858N20151204S0501 trimmed
GIREDUCE: Image tmpfile4365_858N20151204S0502 trimmed
GIREDUCE: Image tmpfile4365_858N20151204S0503 trimmed
GIREDUCE: Image tmpfile4365_858N20151204S0504 trimmed
GIREDUCE: Image tmpfile4365_858N20151204S0505 trimmed


GIREDUCE - Cleaning up -- Sat Apr  9 14:25:35 HST 2016

GIREDUCE -- Sat Apr  9 14:25:35 HST 2016
GIREDUCE exit status:  good.
----------------------------------------------------------------------------
--------------------------------------------------------------------------------
GEMCOMBINE -- Sat Apr  9 14:25:35 HST 2016
 
GEMCOMBINE: input = @tmpcombin4365_855
GEMCOMBINE: output = N20151204S0501_bias.fits
GEMCOMBINE: title = Bias
GEMCOMBINE: bpmfile = 
GEMCOMBINE: sci_ext = SCI
GEMCOMBINE: var_ext = VAR
GEMCOMBINE: dq_ext = DQ
 
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:25: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:25
                Images     Exp
  tmpfile4365_858N20151204S0501.fits[SCI,1]    0.0
  tmpfile4365_858N20151204S0502.fits[SCI,1]    0.0
  tmpfile4365_858N20151204S0503.fits[SCI,1]    0.0
  tmpfile4365_858N20151204S0504.fits[SCI,1]    0.0
  tmpfile4365_858N20151204S0505.fits[SCI,1]    0.0

  Output image = tmpcombout4365_923, ncombine = 5
  Number rejected mask = tmpcombdq4365_924.pl
  Sigma image = tmpcombsig4365_925
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:25: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:25
                Images     Exp
  tmpfile4365_858N20151204S0501.fits[SCI,2]    0.0
  tmpfile4365_858N20151204S0502.fits[SCI,2]    0.0
  tmpfile4365_858N20151204S0503.fits[SCI,2]    0.0
  tmpfile4365_858N20151204S0504.fits[SCI,2]    0.0
  tmpfile4365_858N20151204S0505.fits[SCI,2]    0.0

  Output image = tmpcombout4365_930, ncombine = 5
  Number rejected mask = tmpcombdq4365_931.pl
  Sigma image = tmpcombsig4365_932
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:25: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:25
                Images     Exp
  tmpfile4365_858N20151204S0501.fits[SCI,3]    0.0
  tmpfile4365_858N20151204S0502.fits[SCI,3]    0.0
  tmpfile4365_858N20151204S0503.fits[SCI,3]    0.0
  tmpfile4365_858N20151204S0504.fits[SCI,3]    0.0
  tmpfile4365_858N20151204S0505.fits[SCI,3]    0.0

  Output image = tmpcombout4365_934, ncombine = 5
  Number rejected mask = tmpcombdq4365_935.pl
  Sigma image = tmpcombsig4365_936
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:25: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:25
                Images     Exp
  tmpfile4365_858N20151204S0501.fits[SCI,4]    0.0
  tmpfile4365_858N20151204S0502.fits[SCI,4]    0.0
  tmpfile4365_858N20151204S0503.fits[SCI,4]    0.0
  tmpfile4365_858N20151204S0504.fits[SCI,4]    0.0
  tmpfile4365_858N20151204S0505.fits[SCI,4]    0.0

  Output image = tmpcombout4365_938, ncombine = 5
  Number rejected mask = tmpcombdq4365_939.pl
  Sigma image = tmpcombsig4365_940
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:25: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:25
                Images     Exp
  tmpfile4365_858N20151204S0501.fits[SCI,5]    0.0
  tmpfile4365_858N20151204S0502.fits[SCI,5]    0.0
  tmpfile4365_858N20151204S0503.fits[SCI,5]    0.0
  tmpfile4365_858N20151204S0504.fits[SCI,5]    0.0
  tmpfile4365_858N20151204S0505.fits[SCI,5]    0.0

  Output image = tmpcombout4365_942, ncombine = 5
  Number rejected mask = tmpcombdq4365_943.pl
  Sigma image = tmpcombsig4365_944
WARNING - GEMCOMBINE: five or less images to combine.

Apr  9 14:25: IMCOMBINE
  combine = average, scale = none, zero = none, weight = none
  reject = avsigclip, mclip = yes, nkeep = 1
  lsigma = 3., hsigma = 3.
  blank = 0.
  statsec = Apr  9 14:25
                Images     Exp
  tmpfile4365_858N20151204S0501.fits[SCI,6]    0.0
  tmpfile4365_858N20151204S0502.fits[SCI,6]    0.0
  tmpfile4365_858N20151204S0503.fits[SCI,6]    0.0
  tmpfile4365_858N20151204S0504.fits[SCI,6]    0.0
  tmpfile4365_858N20151204S0505.fits[SCI,6]    0.0

  Output image = tmpcombout4365_946, ncombine = 5
  Number rejected mask = tmpcombdq4365_947.pl
  Sigma image = tmpcombsig4365_948
 
GEMCOMBINE -- Finished: Sat Apr  9 14:25:46 HST 2016
 
GEMCOMBINE -- Exit status: GOOD
--------------------------------------------------------------------------------

GBIAS: Exit status: GOOD.
----------------------------------------------------------------------------
Out[8]:
['N20151204S0501_bias.fits',
 'N20151205S0525_bias.fits',
 'N20151208S0272_bias.fits',
 'N20151210S0430_bias.fits',
 'N20151214S0346_bias.fits']

Visually inspect the processed biases


In [9]:
# WARNING: interactive.  It launches imexam.  Type 'q' to quit.
for biaskey in biases.keys():
    iraf.gdisplay(biases[biaskey]['procbias'], 1, fl_paste='yes')


GDISPLAY - Started: 2016-04-09T14:27:50

GDISPLAY: Temporary image name tmpout4365_951N20151210S0430_bias.fits

GMULTIAMP - Started: 2016-04-09T14:27:50

GMULTIAMP: Tiling image...
               Image tiled.

GMULTIAMP - Finished: 2016-04-09T14:27:52
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY: Starting imexamine - quit with q

GDISPLAY - Finished: 2016-04-09T14:27:58


GDISPLAY - Started: 2016-04-09T14:27:58

GDISPLAY: Temporary image name tmpout4365_985N20151214S0346_bias.fits

GMULTIAMP - Started: 2016-04-09T14:27:58

GMULTIAMP: Tiling image...
               Image tiled.

GMULTIAMP - Finished: 2016-04-09T14:28:01
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY: Starting imexamine - quit with q

GDISPLAY - Finished: 2016-04-09T14:28:02


GDISPLAY - Started: 2016-04-09T14:28:02

GDISPLAY: Temporary image name tmpout4365_1019N20151204S0501_bias.fits

GMULTIAMP - Started: 2016-04-09T14:28:02

GMULTIAMP: Tiling image...
               Image tiled.

GMULTIAMP - Finished: 2016-04-09T14:28:04
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY: Starting imexamine - quit with q

GDISPLAY - Finished: 2016-04-09T14:28:06

STEP 5: Process the twilight flats

This will produce the one twilight flat defined above. That flat will be used later to process the science frames and the fringe frames.


In [11]:
os.chdir(redux_path)

for flatkey in flats.keys():
    
    flatbiassets = zip(flats[flatkey]['frames'], flats[flatkey]['flatbias'])
    for ((root, frames), bias) in flatbiassets:
        flatlist = flatkey+'_tmp.lis'
        iraf.gemlist(root, frames, Stdout=flatlist)
        
        iraf.imdelete('g//@'+flatlist, verify='no')
        iraf.imdelete('rg//@'+flatlist, verify='no')
        
        procbias = os.path.join(procbias_path, bias)
        
        print 'Processing raw twilight flats'
        
        # gpreare and bias correct the raw flats
        iraf.gireduce('@'+flatlist, rawpath=raw_path, \
                  fl_over='yes', fl_trim='yes', \
                  fl_bias='yes', bias=procbias, fl_flat='no', \
                  fl_vardq='yes')
        iraf.delete(flatlist)
        
    
    
#    iraf.imdelete('g//@'+flatkey+'.lis', verify='no')
#    iraf.imdelete('rg//@'+flatkey+'.lis', verify='no')
#    
#    flatlist = flatkey+'.lis'
#    procbias = os.path.join(procbias_path, flats[flatkey]['flatbias'])
#    
#    print 'Processing raw twilight flats'
#    
#    # gprepare and bias correct the raw flats
#    iraf.gireduce('@'+flatlist, rawpath=raw_path, \
#                  fl_over='yes', fl_trim='yes', \
#                  fl_bias='yes', bias=procbias, fl_flat='no', \
#                  fl_vardq='yes')


Processing raw twilight flats
----------------------------------------------------------------------------
GIREDUCE -- Sat Apr  9 14:29:08 HST 2016

----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:29:08 HST 2016

Input list    = N20151204S0446
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0446.fits   Output  gN20151204S0446.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:29:09 HST 2016

Input list    = N20151204S0447
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0447.fits   Output  gN20151204S0447.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:29:10 HST 2016

Input list    = N20151204S0448
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0448.fits   Output  gN20151204S0448.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:29:12 HST 2016

Input list    = N20151204S0449
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0449.fits   Output  gN20151204S0449.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:29:13 HST 2016

Input list    = N20151204S0450
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0450.fits   Output  gN20151204S0450.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:29:14 HST 2016

Input list    = N20151204S0451
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151204S0451.fits   Output  gN20151204S0451.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------

Input files:
gN20151204S0446
gN20151204S0447
gN20151204S0448
gN20151204S0449
gN20151204S0450
gN20151204S0451

Output files:
rgN20151204S0446
rgN20151204S0447
rgN20151204S0448
rgN20151204S0449
rgN20151204S0450
rgN20151204S0451

WARNING - GIREDUCE: Bad Pixel Mask filename is an empty string
                    Only saturated pixels will be flagged
GIREDUCE: Image rgN20151204S0446 overscan subtracted
GIREDUCE: Image rgN20151204S0447 overscan subtracted
GIREDUCE: Image rgN20151204S0448 overscan subtracted
GIREDUCE: Image rgN20151204S0449 overscan subtracted
GIREDUCE: Image rgN20151204S0450 overscan subtracted
GIREDUCE: Image rgN20151204S0451 overscan subtracted
GIREDUCE: Image rgN20151204S0446 trimmed
GIREDUCE: Image rgN20151204S0447 trimmed
GIREDUCE: Image rgN20151204S0448 trimmed
GIREDUCE: Image rgN20151204S0449 trimmed
GIREDUCE: Image rgN20151204S0450 trimmed
GIREDUCE: Image rgN20151204S0451 trimmed

Output image       Bias                 Flat               Dark       Scale 
rgN20151204S0446   /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151204S0501_bias INDEF              INDEF       0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151204S0501_bias from rgN20151204S0446
GIREDUCE: multiplying image rgN20151204S0446 by gain
GIREDUCE: output counts in electrons
Output image       Bias                 Flat               Dark       Scale 
rgN20151204S0447   /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151204S0501_bias INDEF              INDEF       0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151204S0501_bias from rgN20151204S0447
GIREDUCE: multiplying image rgN20151204S0447 by gain
GIREDUCE: output counts in electrons
Output image       Bias                 Flat               Dark       Scale 
rgN20151204S0448   /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151204S0501_bias INDEF              INDEF       0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151204S0501_bias from rgN20151204S0448
GIREDUCE: multiplying image rgN20151204S0448 by gain
GIREDUCE: output counts in electrons
Output image       Bias                 Flat               Dark       Scale 
rgN20151204S0449   /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151204S0501_bias INDEF              INDEF       0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151204S0501_bias from rgN20151204S0449
GIREDUCE: multiplying image rgN20151204S0449 by gain
GIREDUCE: output counts in electrons
Output image       Bias                 Flat               Dark       Scale 
rgN20151204S0450   /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151204S0501_bias INDEF              INDEF       0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151204S0501_bias from rgN20151204S0450
GIREDUCE: multiplying image rgN20151204S0450 by gain
GIREDUCE: output counts in electrons
Output image       Bias                 Flat               Dark       Scale 
rgN20151204S0451   /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151204S0501_bias INDEF              INDEF       0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151204S0501_bias from rgN20151204S0451
GIREDUCE: multiplying image rgN20151204S0451 by gain
GIREDUCE: output counts in electrons

GIREDUCE - Cleaning up -- Sat Apr  9 14:31:28 HST 2016

GIREDUCE -- Sat Apr  9 14:31:28 HST 2016
GIREDUCE exit status:  good.
----------------------------------------------------------------------------
Processing raw twilight flats
----------------------------------------------------------------------------
GIREDUCE -- Sat Apr  9 14:31:29 HST 2016

----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:31:29 HST 2016

Input list    = N20151210S0359
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0359.fits   Output  gN20151210S0359.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:31:30 HST 2016

Input list    = N20151210S0360
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0360.fits   Output  gN20151210S0360.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:31:31 HST 2016

Input list    = N20151210S0361
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151210S0361.fits   Output  gN20151210S0361.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------

Input files:
gN20151210S0359
gN20151210S0360
gN20151210S0361

Output files:
rgN20151210S0359
rgN20151210S0360
rgN20151210S0361

WARNING - GIREDUCE: Bad Pixel Mask filename is an empty string
                    Only saturated pixels will be flagged
GIREDUCE: Image rgN20151210S0359 overscan subtracted
GIREDUCE: Image rgN20151210S0360 overscan subtracted
GIREDUCE: Image rgN20151210S0361 overscan subtracted
GIREDUCE: Image rgN20151210S0359 trimmed
GIREDUCE: Image rgN20151210S0360 trimmed
GIREDUCE: Image rgN20151210S0361 trimmed

Output image       Bias                 Flat               Dark       Scale 
rgN20151210S0359   /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151210S0430_bias INDEF              INDEF       0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151210S0430_bias from rgN20151210S0359
GIREDUCE: multiplying image rgN20151210S0359 by gain
GIREDUCE: output counts in electrons
Output image       Bias                 Flat               Dark       Scale 
rgN20151210S0360   /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151210S0430_bias INDEF              INDEF       0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151210S0430_bias from rgN20151210S0360
GIREDUCE: multiplying image rgN20151210S0360 by gain
GIREDUCE: output counts in electrons
Output image       Bias                 Flat               Dark       Scale 
rgN20151210S0361   /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151210S0430_bias INDEF              INDEF       0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/NGC3310/calib/bias/N20151210S0430_bias from rgN20151210S0361
GIREDUCE: multiplying image rgN20151210S0361 by gain
GIREDUCE: output counts in electrons

GIREDUCE - Cleaning up -- Sat Apr  9 14:32:38 HST 2016

GIREDUCE -- Sat Apr  9 14:32:38 HST 2016
GIREDUCE exit status:  good.
----------------------------------------------------------------------------

In [12]:
os.chdir(redux_path)

for flatkey in flats.keys():
    iraf.imdelete(flats[flatkey]['procflat'], verify='no')
    
    flatlist = flatkey+'.lis'
    procflat = flats[flatkey]['procflat']
    
    print 'Creating twilight flat ', procflat
    
    # Create the processed flat
    iraf.giflat('rg//@'+flatlist, outflat=procflat, fl_vardq='yes')
    
    # Copy the processed flat to the calibration directory
    shutil.copy(procflat, procflat_path)

# The processed flats are:
os.listdir(procflat_path)


Creating twilight flat  N2015120402S0446_flat_i.fits
----------------------------------------------------------------------------
GIFLAT -- Sat Apr  9 14:36:55 HST 2016
 
Flat field images = rg//@flat_20151204.lis
Output flat field image = N2015120402S0446_flat_i.fits
 
WARNING - GIFLAT: using previously gireduce'd image: rgN20151204S0446
WARNING - GIFLAT: using previously gireduce'd image: rgN20151204S0447
WARNING - GIFLAT: using previously gireduce'd image: rgN20151204S0448
WARNING - GIFLAT: using previously gireduce'd image: rgN20151204S0449
WARNING - GIFLAT: using previously gireduce'd image: rgN20151204S0450
WARNING - GIFLAT: using previously gireduce'd image: rgN20151204S0451
WARNING - GIFLAT: using previously gireduce'd image: rgN20151210S0359
WARNING - GIFLAT: using previously gireduce'd image: rgN20151210S0360
WARNING - GIFLAT: using previously gireduce'd image: rgN20151210S0361
Using input files:
rgN20151204S0446
rgN20151204S0447
rgN20151204S0448
rgN20151204S0449
rgN20151204S0450
rgN20151204S0451
rgN20151210S0359
rgN20151210S0360
rgN20151210S0361
Using STATSEC: [50:900,50:2250] from CCD2
Using NORMSEC: [50:900,50:2250] from CCD2
GIFLAT: Normalization factor 110930.8

GIFLAT exit status: good.
----------------------------------------------------------------------------
Out[12]:
['N2015120402S0446_flat_i.fits',
 'N20151205S0018_flat_g.fits',
 'N20151208S0225_flat_r.fits']

Visually inspect the processed twilight flats


In [13]:
for flatkey in flats.keys():
    iraf.gdisplay(flats[flatkey]['procflat'], 1, fl_paste='yes')


GDISPLAY - Started: 2016-04-09T14:37:38

GDISPLAY: Temporary image name tmpout4365_1476N2015120402S0446_flat_i.fits

GMULTIAMP - Started: 2016-04-09T14:37:39

GMULTIAMP: Tiling image...
               Image tiled.

GMULTIAMP - Finished: 2016-04-09T14:37:41
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY: Starting imexamine - quit with q

GDISPLAY - Finished: 2016-04-09T14:37:49

STEP 6: Process the fringe frame

STEP 7: Process the target

In the redder bands, for some CCDs, the fringing can be important. With the EEV CCDs, that is the case for the i-band and the z-band. For this dataset, the CCDs are the newer e2V CCDs and fringing, even in the i-band is minimal such that no correction is necessary.

STEP 7.1: Bias and flat fielding of the target


In [14]:
os.chdir(redux_path)

iraf.imdelete('g//@sci.lis', verify='no')
iraf.imdelete('rg//@sci.lis', verify='no')

iraf.gireduce('@sci.lis', rawpath=raw_path, fl_over='yes', fl_trim='yes', \
              fl_bias='yes', bias=science['scibias'], \
              fl_flat='yes', flat1=science['sciflat'], \
              fl_vardq='yes')


----------------------------------------------------------------------------
GIREDUCE -- Sat Apr  9 14:40:11 HST 2016

----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:40:11 HST 2016

Input list    = N20151214S0239
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0239.fits   Output  gN20151214S0239.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:40:13 HST 2016

Input list    = N20151214S0240
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0240.fits   Output  gN20151214S0240.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:40:14 HST 2016

Input list    = N20151214S0241
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0241.fits   Output  gN20151214S0241.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Sat Apr  9 14:40:15 HST 2016

Input list    = N20151214S0242
Output list   = 
Output prefix = g
Raw path      = /data/workspace/DRWorkshopAAO/NGC3310/raw/
MDF dir       = gmos$data/
Add MDF       = no

Input  /data/workspace/DRWorkshopAAO/NGC3310/raw/N20151214S0242.fits   Output  gN20151214S0242.fits

GPREPARE exit status: good.
----------------------------------------------------------------------------

Input files:
gN20151214S0239
gN20151214S0240
gN20151214S0241
gN20151214S0242

Output files:
rgN20151214S0239
rgN20151214S0240
rgN20151214S0241
rgN20151214S0242

WARNING - GIREDUCE: Bad Pixel Mask filename is an empty string
                    Only saturated pixels will be flagged
GIREDUCE: Image rgN20151214S0239 overscan subtracted
GIREDUCE: Image rgN20151214S0240 overscan subtracted
GIREDUCE: Image rgN20151214S0241 overscan subtracted
GIREDUCE: Image rgN20151214S0242 overscan subtracted
GIREDUCE: Image rgN20151214S0239 trimmed
GIREDUCE: Image rgN20151214S0240 trimmed
GIREDUCE: Image rgN20151214S0241 trimmed
GIREDUCE: Image rgN20151214S0242 trimmed

Output image       Bias                 Flat               Dark       Scale 
rgN20151214S0239   N20151214S0346_bias  N2015120402S0446_flat_i INDEF       0.00
GIREDUCE: Subtracted bias N20151214S0346_bias from rgN20151214S0239
GIREDUCE: Divided image rgN20151214S0239 by flatfield N2015120402S0446_flat_i
GIREDUCE: multiplying image rgN20151214S0239 by gain
GIREDUCE: output counts in electrons
Output image       Bias                 Flat               Dark       Scale 
rgN20151214S0240   N20151214S0346_bias  N2015120402S0446_flat_i INDEF       0.00
GIREDUCE: Subtracted bias N20151214S0346_bias from rgN20151214S0240
GIREDUCE: Divided image rgN20151214S0240 by flatfield N2015120402S0446_flat_i
GIREDUCE: multiplying image rgN20151214S0240 by gain
GIREDUCE: output counts in electrons
Output image       Bias                 Flat               Dark       Scale 
rgN20151214S0241   N20151214S0346_bias  N2015120402S0446_flat_i INDEF       0.00
GIREDUCE: Subtracted bias N20151214S0346_bias from rgN20151214S0241
GIREDUCE: Divided image rgN20151214S0241 by flatfield N2015120402S0446_flat_i
GIREDUCE: multiplying image rgN20151214S0241 by gain
GIREDUCE: output counts in electrons
Output image       Bias                 Flat               Dark       Scale 
rgN20151214S0242   N20151214S0346_bias  N2015120402S0446_flat_i INDEF       0.00
GIREDUCE: Subtracted bias N20151214S0346_bias from rgN20151214S0242
GIREDUCE: Divided image rgN20151214S0242 by flatfield N2015120402S0446_flat_i
GIREDUCE: multiplying image rgN20151214S0242 by gain
GIREDUCE: output counts in electrons

GIREDUCE - Cleaning up -- Sat Apr  9 14:41:52 HST 2016

GIREDUCE -- Sat Apr  9 14:41:52 HST 2016
GIREDUCE exit status:  good.
----------------------------------------------------------------------------

Visually inspect the reduce science frames


In [15]:
sciframes = open('sci.lis', 'r')
for line in sciframes:
    frame = line.strip()
    iraf.gdisplay('rg'+frame, 1, fl_paste='yes')
sciframes.close()


GDISPLAY - Started: 2016-04-09T14:42:00

GDISPLAY: Temporary image name tmpout4365_1604rgN20151214S0239.fits

GMULTIAMP - Started: 2016-04-09T14:42:00

GMULTIAMP: Tiling image...
               Image tiled.

GMULTIAMP - Finished: 2016-04-09T14:42:02
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY: Starting imexamine - quit with q

GDISPLAY - Finished: 2016-04-09T14:42:19


GDISPLAY - Started: 2016-04-09T14:42:19

GDISPLAY: Temporary image name tmpout4365_1638rgN20151214S0240.fits

GMULTIAMP - Started: 2016-04-09T14:42:19

GMULTIAMP: Tiling image...
               Image tiled.

GMULTIAMP - Finished: 2016-04-09T14:42:21
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY: Starting imexamine - quit with q

GDISPLAY - Finished: 2016-04-09T14:42:23


GDISPLAY - Started: 2016-04-09T14:42:23

GDISPLAY: Temporary image name tmpout4365_1672rgN20151214S0241.fits

GMULTIAMP - Started: 2016-04-09T14:42:23

GMULTIAMP: Tiling image...
               Image tiled.

GMULTIAMP - Finished: 2016-04-09T14:42:25
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY: Starting imexamine - quit with q

GDISPLAY - Finished: 2016-04-09T14:42:26


GDISPLAY - Started: 2016-04-09T14:42:26

GDISPLAY: Temporary image name tmpout4365_1706rgN20151214S0242.fits

GMULTIAMP - Started: 2016-04-09T14:42:27

GMULTIAMP: Tiling image...
               Image tiled.

GMULTIAMP - Finished: 2016-04-09T14:42:29
GMULTIAMP - Exit status: GOOD

GDISPLAY: Displaying image...

GDISPLAY: Starting imexamine - quit with q

GDISPLAY - Finished: 2016-04-09T14:42:30

STEP 7.2: Fringe removal

STEP 8: Mosaic the reduced target frames


In [16]:
iraf.imdelete('mrg//@sci.lis', verify='no')

iraf.gmosaic('rg//@sci.lis', fl_vardq='yes')


GMOSAIC 
GMOSAIC -------------------------------------------------------------------
GMOSAIC Log opened at [Sat 14:42:37 09-Apr-2016]
GMOSAIC 
GMOSAIC --------------------
GMOSAIC Input Parameters:
GMOSAIC      inimages       = rg//@sci.lis
GMOSAIC      outimages      = 
GMOSAIC      outpref        = m
GMOSAIC      fl_paste       = no
GMOSAIC      fl_vardq       = yes
GMOSAIC      fl_fixpix      = no
GMOSAIC      fl_clean       = yes
GMOSAIC      fl_fulldq      = no
GMOSAIC      bitflags       = all
GMOSAIC      dqthresh       = 0.1
GMOSAIC      geointer       = linear
GMOSAIC      gap            = default
GMOSAIC      bpmfile        = gmos$data/chipgaps.dat
GMOSAIC      statsec        = default
GMOSAIC      obsmode        = IMAGE
GMOSAIC      sci_ext        = SCI
GMOSAIC      var_ext        = VAR
GMOSAIC      dq_ext         = DQ
GMOSAIC      key_detsec     = DETSEC
GMOSAIC      key_ccdsec     = CCDSEC
GMOSAIC      key_datsec     = DATASEC
GMOSAIC      key_ccdsum     = CCDSUM
GMOSAIC      key_obsmode    = OBSMODE
GMOSAIC      logfile        = 
GMOSAIC      verbose        = yes
GMOSAIC      fl_real        = no
GMOSAIC --------------------
GMOSAIC Input: rgN20151214S0239  Output: mrgN20151214S0239.fits
GMOSAIC imcopy  rgN20151214S0239[SCI,2][1:512,1:2304]  tmptile4365_1755[1,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[SCI,1][1:512,1:2304]  tmptile4365_1755[1,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[VAR,2][1:512,1:2304]  tmptile4365_1755[2,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[VAR,1][1:512,1:2304]  tmptile4365_1755[2,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[DQ,2][1:512,1:2304]  tmptile4365_1755[3,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[DQ,1][1:512,1:2304]  tmptile4365_1755[3,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[SCI,4][1:512,1:2304]  tmptile4365_1756[1,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[SCI,3][1:512,1:2304]  tmptile4365_1756[1,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[VAR,4][1:512,1:2304]  tmptile4365_1756[2,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[VAR,3][1:512,1:2304]  tmptile4365_1756[2,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[DQ,4][1:512,1:2304]  tmptile4365_1756[3,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[DQ,3][1:512,1:2304]  tmptile4365_1756[3,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[SCI,5][1:512,1:2304]  tmptile4365_1757[1,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[SCI,6][1:512,1:2304]  tmptile4365_1757[1,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[VAR,5][1:512,1:2304]  tmptile4365_1757[2,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[VAR,6][1:512,1:2304]  tmptile4365_1757[2,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[DQ,5][1:512,1:2304]  tmptile4365_1757[3,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0239[DQ,6][1:512,1:2304]  tmptile4365_1757[3,append][513:1024,1:2304]  INDEF  0.
GMOSAIC Setting rotation Xrot[1]=-0.009000, Yrot[1]=-0.009000
GMOSAIC Setting rotation Xrot[2]= 0.000000, Yrot[2]= 0.000000
GMOSAIC Setting rotation Xrot[3]=-0.003000, Yrot[3]=-0.003000
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1751[SCI,1][1:1024,1:2304] to image tmpchipsci4365_1758
GMOSAIC     xshift: 2.03 yshift: 0.29 xmag: 1.00 ymag: 1.00 xrot: 359.99 yrot: 359.99
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1751[VAR,1][1:1024,1:2304] to image tmpchipvar4365_1759
GMOSAIC     xshift: 2.03 yshift: 0.29 xmag: 1.00 ymag: 1.00 xrot: 359.99 yrot: 359.99
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1751[DQ,1][1:1024,1:2304] to image tmpdq_transimg4365_1747.fits
GMOSAIC     xshift: 2.03 yshift: 0.29 xmag: 1.00 ymag: 1.00 xrot: 359.99 yrot: 359.99
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1751[SCI,2][1:1024,1:2304] to image tmpchipsci4365_1761
GMOSAIC     xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1751[VAR,2][1:1024,1:2304] to image tmpchipvar4365_1762
GMOSAIC     xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1751[SCI,3][1:1024,1:2304] to image tmpchipsci4365_1764
GMOSAIC     xshift: -1.84 yshift: -1.05 xmag: 1.00 ymag: 1.00 xrot: 360.00 yrot: 360.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1751[VAR,3][1:1024,1:2304] to image tmpchipvar4365_1765
GMOSAIC     xshift: -1.84 yshift: -1.05 xmag: 1.00 ymag: 1.00 xrot: 360.00 yrot: 360.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1751[DQ,3][1:1024,1:2304] to image tmpdq_transimg4365_1747.fits
GMOSAIC     xshift: -1.84 yshift: -1.05 xmag: 1.00 ymag: 1.00 xrot: 360.00 yrot: 360.00
GMOSAIC imcopy  tmpchipsci4365_1758  mrgN20151214S0239.fits[1,append][1:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipsci4365_1761  mrgN20151214S0239.fits[1,append][1043:2066,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipsci4365_1764  mrgN20151214S0239.fits[1,append][2085:3108,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipvar4365_1759  mrgN20151214S0239.fits[2,append][1:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipvar4365_1762  mrgN20151214S0239.fits[2,append][1043:2066,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipvar4365_1765  mrgN20151214S0239.fits[2,append][2085:3108,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipdq4365_1760  mrgN20151214S0239.fits[3,append][1:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipdq4365_1763  mrgN20151214S0239.fits[3,append][1043:2066,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipdq4365_1766  mrgN20151214S0239.fits[3,append][2085:3108,1:2304]  INDEF  0.
gmos$data/gmos-n_bpm_e2v_22_6amp_v1_mosaic.pl
GMOSAIC Cleaning areas outside imaging field, cleaning value = 3037.511
GMOSAIC Input: rgN20151214S0240  Output: mrgN20151214S0240.fits
GMOSAIC imcopy  rgN20151214S0240[SCI,2][1:512,1:2304]  tmptile4365_1771[1,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[SCI,1][1:512,1:2304]  tmptile4365_1771[1,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[VAR,2][1:512,1:2304]  tmptile4365_1771[2,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[VAR,1][1:512,1:2304]  tmptile4365_1771[2,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[DQ,2][1:512,1:2304]  tmptile4365_1771[3,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[DQ,1][1:512,1:2304]  tmptile4365_1771[3,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[SCI,4][1:512,1:2304]  tmptile4365_1772[1,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[SCI,3][1:512,1:2304]  tmptile4365_1772[1,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[VAR,4][1:512,1:2304]  tmptile4365_1772[2,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[VAR,3][1:512,1:2304]  tmptile4365_1772[2,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[DQ,4][1:512,1:2304]  tmptile4365_1772[3,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[DQ,3][1:512,1:2304]  tmptile4365_1772[3,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[SCI,5][1:512,1:2304]  tmptile4365_1773[1,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[SCI,6][1:512,1:2304]  tmptile4365_1773[1,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[VAR,5][1:512,1:2304]  tmptile4365_1773[2,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[VAR,6][1:512,1:2304]  tmptile4365_1773[2,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[DQ,5][1:512,1:2304]  tmptile4365_1773[3,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0240[DQ,6][1:512,1:2304]  tmptile4365_1773[3,append][513:1024,1:2304]  INDEF  0.
GMOSAIC Setting rotation Xrot[1]=-0.009000, Yrot[1]=-0.009000
GMOSAIC Setting rotation Xrot[2]= 0.000000, Yrot[2]= 0.000000
GMOSAIC Setting rotation Xrot[3]=-0.003000, Yrot[3]=-0.003000
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1767[SCI,1][1:1024,1:2304] to image tmpchipsci4365_1774
GMOSAIC     xshift: 2.03 yshift: 0.29 xmag: 1.00 ymag: 1.00 xrot: 359.99 yrot: 359.99
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1767[VAR,1][1:1024,1:2304] to image tmpchipvar4365_1775
GMOSAIC     xshift: 2.03 yshift: 0.29 xmag: 1.00 ymag: 1.00 xrot: 359.99 yrot: 359.99
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1767[DQ,1][1:1024,1:2304] to image tmpdq_transimg4365_1747.fits
GMOSAIC     xshift: 2.03 yshift: 0.29 xmag: 1.00 ymag: 1.00 xrot: 359.99 yrot: 359.99
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1767[SCI,2][1:1024,1:2304] to image tmpchipsci4365_1777
GMOSAIC     xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1767[VAR,2][1:1024,1:2304] to image tmpchipvar4365_1778
GMOSAIC     xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1767[SCI,3][1:1024,1:2304] to image tmpchipsci4365_1780
GMOSAIC     xshift: -1.84 yshift: -1.05 xmag: 1.00 ymag: 1.00 xrot: 360.00 yrot: 360.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1767[VAR,3][1:1024,1:2304] to image tmpchipvar4365_1781
GMOSAIC     xshift: -1.84 yshift: -1.05 xmag: 1.00 ymag: 1.00 xrot: 360.00 yrot: 360.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1767[DQ,3][1:1024,1:2304] to image tmpdq_transimg4365_1747.fits
GMOSAIC     xshift: -1.84 yshift: -1.05 xmag: 1.00 ymag: 1.00 xrot: 360.00 yrot: 360.00
GMOSAIC imcopy  tmpchipsci4365_1774  mrgN20151214S0240.fits[1,append][1:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipsci4365_1777  mrgN20151214S0240.fits[1,append][1043:2066,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipsci4365_1780  mrgN20151214S0240.fits[1,append][2085:3108,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipvar4365_1775  mrgN20151214S0240.fits[2,append][1:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipvar4365_1778  mrgN20151214S0240.fits[2,append][1043:2066,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipvar4365_1781  mrgN20151214S0240.fits[2,append][2085:3108,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipdq4365_1776  mrgN20151214S0240.fits[3,append][1:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipdq4365_1779  mrgN20151214S0240.fits[3,append][1043:2066,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipdq4365_1782  mrgN20151214S0240.fits[3,append][2085:3108,1:2304]  INDEF  0.
gmos$data/gmos-n_bpm_e2v_22_6amp_v1_mosaic.pl
GMOSAIC Cleaning areas outside imaging field, cleaning value = 3067.38
GMOSAIC Input: rgN20151214S0241  Output: mrgN20151214S0241.fits
GMOSAIC imcopy  rgN20151214S0241[SCI,2][1:512,1:2304]  tmptile4365_1787[1,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[SCI,1][1:512,1:2304]  tmptile4365_1787[1,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[VAR,2][1:512,1:2304]  tmptile4365_1787[2,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[VAR,1][1:512,1:2304]  tmptile4365_1787[2,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[DQ,2][1:512,1:2304]  tmptile4365_1787[3,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[DQ,1][1:512,1:2304]  tmptile4365_1787[3,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[SCI,4][1:512,1:2304]  tmptile4365_1788[1,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[SCI,3][1:512,1:2304]  tmptile4365_1788[1,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[VAR,4][1:512,1:2304]  tmptile4365_1788[2,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[VAR,3][1:512,1:2304]  tmptile4365_1788[2,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[DQ,4][1:512,1:2304]  tmptile4365_1788[3,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[DQ,3][1:512,1:2304]  tmptile4365_1788[3,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[SCI,5][1:512,1:2304]  tmptile4365_1789[1,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[SCI,6][1:512,1:2304]  tmptile4365_1789[1,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[VAR,5][1:512,1:2304]  tmptile4365_1789[2,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[VAR,6][1:512,1:2304]  tmptile4365_1789[2,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[DQ,5][1:512,1:2304]  tmptile4365_1789[3,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0241[DQ,6][1:512,1:2304]  tmptile4365_1789[3,append][513:1024,1:2304]  INDEF  0.
GMOSAIC Setting rotation Xrot[1]=-0.009000, Yrot[1]=-0.009000
GMOSAIC Setting rotation Xrot[2]= 0.000000, Yrot[2]= 0.000000
GMOSAIC Setting rotation Xrot[3]=-0.003000, Yrot[3]=-0.003000
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1783[SCI,1][1:1024,1:2304] to image tmpchipsci4365_1790
GMOSAIC     xshift: 2.03 yshift: 0.29 xmag: 1.00 ymag: 1.00 xrot: 359.99 yrot: 359.99
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1783[VAR,1][1:1024,1:2304] to image tmpchipvar4365_1791
GMOSAIC     xshift: 2.03 yshift: 0.29 xmag: 1.00 ymag: 1.00 xrot: 359.99 yrot: 359.99
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1783[DQ,1][1:1024,1:2304] to image tmpdq_transimg4365_1747.fits
GMOSAIC     xshift: 2.03 yshift: 0.29 xmag: 1.00 ymag: 1.00 xrot: 359.99 yrot: 359.99
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1783[SCI,2][1:1024,1:2304] to image tmpchipsci4365_1793
GMOSAIC     xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1783[VAR,2][1:1024,1:2304] to image tmpchipvar4365_1794
GMOSAIC     xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1783[SCI,3][1:1024,1:2304] to image tmpchipsci4365_1796
GMOSAIC     xshift: -1.84 yshift: -1.05 xmag: 1.00 ymag: 1.00 xrot: 360.00 yrot: 360.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1783[VAR,3][1:1024,1:2304] to image tmpchipvar4365_1797
GMOSAIC     xshift: -1.84 yshift: -1.05 xmag: 1.00 ymag: 1.00 xrot: 360.00 yrot: 360.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1783[DQ,3][1:1024,1:2304] to image tmpdq_transimg4365_1747.fits
GMOSAIC     xshift: -1.84 yshift: -1.05 xmag: 1.00 ymag: 1.00 xrot: 360.00 yrot: 360.00
GMOSAIC imcopy  tmpchipsci4365_1790  mrgN20151214S0241.fits[1,append][1:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipsci4365_1793  mrgN20151214S0241.fits[1,append][1043:2066,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipsci4365_1796  mrgN20151214S0241.fits[1,append][2085:3108,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipvar4365_1791  mrgN20151214S0241.fits[2,append][1:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipvar4365_1794  mrgN20151214S0241.fits[2,append][1043:2066,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipvar4365_1797  mrgN20151214S0241.fits[2,append][2085:3108,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipdq4365_1792  mrgN20151214S0241.fits[3,append][1:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipdq4365_1795  mrgN20151214S0241.fits[3,append][1043:2066,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipdq4365_1798  mrgN20151214S0241.fits[3,append][2085:3108,1:2304]  INDEF  0.
gmos$data/gmos-n_bpm_e2v_22_6amp_v1_mosaic.pl
GMOSAIC Cleaning areas outside imaging field, cleaning value = 3083.523
GMOSAIC Input: rgN20151214S0242  Output: mrgN20151214S0242.fits
GMOSAIC imcopy  rgN20151214S0242[SCI,2][1:512,1:2304]  tmptile4365_1803[1,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[SCI,1][1:512,1:2304]  tmptile4365_1803[1,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[VAR,2][1:512,1:2304]  tmptile4365_1803[2,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[VAR,1][1:512,1:2304]  tmptile4365_1803[2,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[DQ,2][1:512,1:2304]  tmptile4365_1803[3,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[DQ,1][1:512,1:2304]  tmptile4365_1803[3,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[SCI,4][1:512,1:2304]  tmptile4365_1804[1,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[SCI,3][1:512,1:2304]  tmptile4365_1804[1,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[VAR,4][1:512,1:2304]  tmptile4365_1804[2,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[VAR,3][1:512,1:2304]  tmptile4365_1804[2,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[DQ,4][1:512,1:2304]  tmptile4365_1804[3,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[DQ,3][1:512,1:2304]  tmptile4365_1804[3,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[SCI,5][1:512,1:2304]  tmptile4365_1805[1,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[SCI,6][1:512,1:2304]  tmptile4365_1805[1,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[VAR,5][1:512,1:2304]  tmptile4365_1805[2,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[VAR,6][1:512,1:2304]  tmptile4365_1805[2,append][513:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[DQ,5][1:512,1:2304]  tmptile4365_1805[3,append][1:512,1:2304]  INDEF  0.
GMOSAIC imcopy  rgN20151214S0242[DQ,6][1:512,1:2304]  tmptile4365_1805[3,append][513:1024,1:2304]  INDEF  0.
GMOSAIC Setting rotation Xrot[1]=-0.009000, Yrot[1]=-0.009000
GMOSAIC Setting rotation Xrot[2]= 0.000000, Yrot[2]= 0.000000
GMOSAIC Setting rotation Xrot[3]=-0.003000, Yrot[3]=-0.003000
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1799[SCI,1][1:1024,1:2304] to image tmpchipsci4365_1806
GMOSAIC     xshift: 2.03 yshift: 0.29 xmag: 1.00 ymag: 1.00 xrot: 359.99 yrot: 359.99
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1799[VAR,1][1:1024,1:2304] to image tmpchipvar4365_1807
GMOSAIC     xshift: 2.03 yshift: 0.29 xmag: 1.00 ymag: 1.00 xrot: 359.99 yrot: 359.99
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1799[DQ,1][1:1024,1:2304] to image tmpdq_transimg4365_1747.fits
GMOSAIC     xshift: 2.03 yshift: 0.29 xmag: 1.00 ymag: 1.00 xrot: 359.99 yrot: 359.99
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1799[SCI,2][1:1024,1:2304] to image tmpchipsci4365_1809
GMOSAIC     xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1799[VAR,2][1:1024,1:2304] to image tmpchipvar4365_1810
GMOSAIC     xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1799[SCI,3][1:1024,1:2304] to image tmpchipsci4365_1812
GMOSAIC     xshift: -1.84 yshift: -1.05 xmag: 1.00 ymag: 1.00 xrot: 360.00 yrot: 360.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1799[VAR,3][1:1024,1:2304] to image tmpchipvar4365_1813
GMOSAIC     xshift: -1.84 yshift: -1.05 xmag: 1.00 ymag: 1.00 xrot: 360.00 yrot: 360.00
GMOSAIC 
GMOSAIC Transforming image tmpimg4365_1799[DQ,3][1:1024,1:2304] to image tmpdq_transimg4365_1747.fits
GMOSAIC     xshift: -1.84 yshift: -1.05 xmag: 1.00 ymag: 1.00 xrot: 360.00 yrot: 360.00
GMOSAIC imcopy  tmpchipsci4365_1806  mrgN20151214S0242.fits[1,append][1:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipsci4365_1809  mrgN20151214S0242.fits[1,append][1043:2066,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipsci4365_1812  mrgN20151214S0242.fits[1,append][2085:3108,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipvar4365_1807  mrgN20151214S0242.fits[2,append][1:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipvar4365_1810  mrgN20151214S0242.fits[2,append][1043:2066,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipvar4365_1813  mrgN20151214S0242.fits[2,append][2085:3108,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipdq4365_1808  mrgN20151214S0242.fits[3,append][1:1024,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipdq4365_1811  mrgN20151214S0242.fits[3,append][1043:2066,1:2304]  INDEF  0.
GMOSAIC imcopy  tmpchipdq4365_1814  mrgN20151214S0242.fits[3,append][2085:3108,1:2304]  INDEF  0.
gmos$data/gmos-n_bpm_e2v_22_6amp_v1_mosaic.pl
GMOSAIC Cleaning areas outside imaging field, cleaning value = 3058.761
GMOSAIC All 4 images successfully mosaiced.
GMOSAIC --------------------
GMOSAIC Exit status: SUCCESS
GMOSAIC Log closed at [Sat 14:43:16 09-Apr-2016]
GMOSAIC -------------------------------------------------------------------
GMOSAIC 

Visually inspect the mosaiced frames


In [17]:
sciframes = open('sci.lis', 'r')
for line in sciframes:
    frame = line.strip()
    iraf.gdisplay('mrg'+frame, 1, z2=20000, fl_paste='yes')
sciframes.close()


GDISPLAY - Started: 2016-04-09T14:43:41

WARNING - GDISPLAY: only one science extension found, turning paste off

WARNING - GDISPLAY: Input image has been processed by GMOSAIC.
                    Switching off the following flags:
                          ret_roi=no

GDISPLAY: Displaying image...

GDISPLAY: Starting imexamine - quit with q

GDISPLAY - Finished: 2016-04-09T14:43:50


GDISPLAY - Started: 2016-04-09T14:43:51

WARNING - GDISPLAY: only one science extension found, turning paste off

WARNING - GDISPLAY: Input image has been processed by GMOSAIC.
                    Switching off the following flags:
                          ret_roi=no

GDISPLAY: Displaying image...

GDISPLAY: Starting imexamine - quit with q

GDISPLAY - Finished: 2016-04-09T14:43:52


GDISPLAY - Started: 2016-04-09T14:43:52

WARNING - GDISPLAY: only one science extension found, turning paste off

WARNING - GDISPLAY: Input image has been processed by GMOSAIC.
                    Switching off the following flags:
                          ret_roi=no

GDISPLAY: Displaying image...

GDISPLAY: Starting imexamine - quit with q

GDISPLAY - Finished: 2016-04-09T14:43:54


GDISPLAY - Started: 2016-04-09T14:43:54

WARNING - GDISPLAY: only one science extension found, turning paste off

WARNING - GDISPLAY: Input image has been processed by GMOSAIC.
                    Switching off the following flags:
                          ret_roi=no

GDISPLAY: Displaying image...

GDISPLAY: Starting imexamine - quit with q

GDISPLAY - Finished: 2016-04-09T14:43:55

STEP 9: Align and stack the mosaic target frames


In [18]:
iraf.imdelete(science['procsci'], verify='no')

iraf.imcoadd('mrg//@sci.lis', outimage=science['procsci'], logfile=logfile)


-----------------------------------------------------------------------------
IMCOADD -- Sat Apr  9 14:44:01 HST 2016

Images (and masks) in list
mrgN20151214S0239
mrgN20151214S0240
mrgN20151214S0241
mrgN20151214S0242
Alignment method: wcs
Making individual masks from DQ
Fixpix input images using DQ
Instrument + camera      : GMOS-N  
Telescope ISS port number: 5
Bad pixel file gmos$data/gmos-n_bpm_e2v_22_6amp_v1_mosaic.pl
Read noise:   3.20 e-   Gain:   1.00 e-/ADU
Co-added cleaned output image: ngc3310_i
Pixel scale: 0.1458
Statistics section: [100:3008,100:2204]
Finding objects in mrgN20151214S0239
Sampling for geotran  : 20  20
Block size for geotran: 2048  2048
Centering objects in mrgN20151214S0239
Number of objects successfully centered in reference image 310
Alignment method: wcs

Coordinate list: imxymatch.1  Transform: mrgN20151214S0240_trn
    Results file: 
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 0.001294235  0.001326204
Coordinate mapping parameters
    Mean Xref and Yref: 1554.5  1152.5
    Mean Xin and Yin: 1527.023  1111.346
    X and Y shift: -27.48515  -41.1475  (xin  yin)
    X and Y scale: 1.000001  1.000001  (xin / xref  yin / yref)
    X and Y axis rotation: 0.00030  0.00030  (degrees  degrees)
Centering objects in mrgN20151214S0240

Xmax =  3108   Ymax = 2304
Entering geomap
Fitting geometry rscale
Iterating a maximum of 5 times

Coordinate list: mrgN20151214S0240_trn  Transform: mrgN20151214S0240_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 11.26623  10.84043
Coordinate mapping parameters
    Mean Xref and Yref: 1610.564  1190.636
    Mean Xin and Yin: 1583.295  1150.927
    X and Y shift: -23.43967  -36.84631  (xin  yin)
    X and Y scale: 0.9978665  0.9978665  (xin / xref  yin / yref)
    X and Y axis rotation: 359.99283  359.99283  (degrees  degrees)

Coordinate list: mrgN20151214S0240_trn  Transform: mrgN20151214S0240_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 4.548396  3.923584
Coordinate mapping parameters
    Mean Xref and Yref: 1618.191  1208.344
    Mean Xin and Yin: 1590.793  1168.508
    X and Y shift: -26.24071  -39.37842  (xin  yin)
    X and Y scale: 0.9993838  0.9993837  (xin / xref  yin / yref)
    X and Y axis rotation: 359.98468  359.98468  (degrees  degrees)

Coordinate list: mrgN20151214S0240_trn  Transform: mrgN20151214S0240_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 0.2454656  0.2507814
Coordinate mapping parameters
    Mean Xref and Yref: 1605.257  1241.333
    Mean Xin and Yin: 1577.828  1200.966
    X and Y shift: -27.32797  -40.26001  (xin  yin)
    X and Y scale: 0.9999408  0.9999408  (xin / xref  yin / yref)
    X and Y axis rotation: 359.99945  359.99945  (degrees  degrees)

Coordinate list: mrgN20151214S0240_trn  Transform: mrgN20151214S0240_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 0.1781242  0.206538
Coordinate mapping parameters
    Mean Xref and Yref: 1639.448  1208.81
    Mean Xin and Yin: 1612.014  1168.477
    X and Y shift: -27.35039  -40.25781  (xin  yin)
    X and Y scale: 0.9999438  0.9999438  (xin / xref  yin / yref)
    X and Y axis rotation: 0.00035  0.00035  (degrees  degrees)

Coordinate list: imxymatch.1  Transform: mrgN20151214S0241_trn
    Results file: 
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 0.001205312  0.001383584
Coordinate mapping parameters
    Mean Xref and Yref: 1554.5  1152.5
    Mean Xin and Yin: 1609.385  1152.457
    X and Y shift: 54.89256  -0.05618692  (xin  yin)
    X and Y scale: 1.000001  1.000001  (xin / xref  yin / yref)
    X and Y axis rotation: 359.99954  359.99954  (degrees  degrees)
Centering objects in mrgN20151214S0241

Xmax =  3108   Ymax = 2304
Entering geomap
Fitting geometry rscale
Iterating a maximum of 5 times

Coordinate list: mrgN20151214S0241_trn  Transform: mrgN20151214S0241_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 10.76837  11.90341
Coordinate mapping parameters
    Mean Xref and Yref: 1617.805  1185.021
    Mean Xin and Yin: 1643.918  1144.155
    X and Y shift: 30.09101  -37.84436  (xin  yin)
    X and Y scale: 0.9981529  0.9981529  (xin / xref  yin / yref)
    X and Y axis rotation: 359.98904  359.98904  (degrees  degrees)

Coordinate list: mrgN20151214S0241_trn  Transform: mrgN20151214S0241_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 5.118413  4.65624
Coordinate mapping parameters
    Mean Xref and Yref: 1625.856  1197.205
    Mean Xin and Yin: 1653.336  1157.12
    X and Y shift: 28.48309  -40.08447  (xin  yin)
    X and Y scale: 0.999468  0.999468  (xin / xref  yin / yref)
    X and Y axis rotation: 359.97498  359.97498  (degrees  degrees)

Coordinate list: mrgN20151214S0241_trn  Transform: mrgN20151214S0241_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 0.2637238  0.3871498
Coordinate mapping parameters
    Mean Xref and Yref: 1609.351  1243.508
    Mean Xin and Yin: 1636.551  1202.455
    X and Y shift: 27.29078  -41.26697  (xin  yin)
    X and Y scale: 0.999995  0.999995  (xin / xref  yin / yref)
    X and Y axis rotation: 359.99725  359.99725  (degrees  degrees)

Coordinate list: mrgN20151214S0241_trn  Transform: mrgN20151214S0241_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 0.1458438  0.1975279
Coordinate mapping parameters
    Mean Xref and Yref: 1613.313  1230.582
    Mean Xin and Yin: 1640.554  1189.387
    X and Y shift: 27.31236  -41.19324  (xin  yin)
    X and Y scale: 0.9999828  0.9999828  (xin / xref  yin / yref)
    X and Y axis rotation: 359.99832  359.99832  (degrees  degrees)

Coordinate list: mrgN20151214S0241_trn  Transform: mrgN20151214S0241_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 0.142934  0.1598223
Coordinate mapping parameters
    Mean Xref and Yref: 1618.731  1257.578
    Mean Xin and Yin: 1645.982  1216.422
    X and Y shift: 27.34186  -41.18835  (xin  yin)
    X and Y scale: 0.9999757  0.9999757  (xin / xref  yin / yref)
    X and Y axis rotation: 359.99765  359.99765  (degrees  degrees)

Coordinate list: imxymatch.1  Transform: mrgN20151214S0242_trn
    Results file: 
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 0.001294721  0.00128864
Coordinate mapping parameters
    Mean Xref and Yref: 1554.5  1152.5
    Mean Xin and Yin: 1581.979  1193.654
    X and Y shift: 27.46739  41.16833  (xin  yin)
    X and Y scale: 1.000001  1.000001  (xin / xref  yin / yref)
    X and Y axis rotation: 0.00055  0.00055  (degrees  degrees)
Centering objects in mrgN20151214S0242

Xmax =  3108   Ymax = 2304
Entering geomap
Fitting geometry rscale
Iterating a maximum of 5 times

Coordinate list: mrgN20151214S0242_trn  Transform: mrgN20151214S0242_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 11.07054  10.96255
Coordinate mapping parameters
    Mean Xref and Yref: 1545.081  1123.439
    Mean Xin and Yin: 1599.878  1125.428
    X and Y shift: 60.35363  1.635742  (xin  yin)
    X and Y scale: 0.9980758  0.9980758  (xin / xref  yin / yref)
    X and Y axis rotation: 359.90884  359.90884  (degrees  degrees)

Coordinate list: mrgN20151214S0242_trn  Transform: mrgN20151214S0242_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 5.209299  4.387035
Coordinate mapping parameters
    Mean Xref and Yref: 1563.415  1129.397
    Mean Xin and Yin: 1619.12  1130.919
    X and Y shift: 56.0129  0.9387207  (xin  yin)
    X and Y scale: 0.9994962  0.9994963  (xin / xref  yin / yref)
    X and Y axis rotation: 359.97021  359.97021  (degrees  degrees)

Coordinate list: mrgN20151214S0242_trn  Transform: mrgN20151214S0242_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 0.7888601  0.3982207
Coordinate mapping parameters
    Mean Xref and Yref: 1535.687  1165.696
    Mean Xin and Yin: 1590.559  1166.144
    X and Y shift: 55.20993  0.4134521  (xin  yin)
    X and Y scale: 0.9998962  0.9998962  (xin / xref  yin / yref)
    X and Y axis rotation: 359.98981  359.98981  (degrees  degrees)

Coordinate list: mrgN20151214S0242_trn  Transform: mrgN20151214S0242_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 0.208289  0.2289933
Coordinate mapping parameters
    Mean Xref and Yref: 1614.449  1206.717
    Mean Xin and Yin: 1669.226  1207.285
    X and Y shift: 54.9104  0.6104736  (xin  yin)
    X and Y scale: 0.9999383  0.9999383  (xin / xref  yin / yref)
    X and Y axis rotation: 359.99820  359.99820  (degrees  degrees)

Coordinate list: mrgN20151214S0242_trn  Transform: mrgN20151214S0242_trn
    Results file: tmpres4365_1836
Coordinate mapping status
    X fit ok.  Y fit ok.
    Xin and Yin fit rms: 0.1786358  0.1967361
Coordinate mapping parameters
    Mean Xref and Yref: 1627.382  1217.822
    Mean Xin and Yin: 1682.166  1218.414
    X and Y shift: 54.91008  0.6824951  (xin  yin)
    X and Y scale: 0.9999247  0.9999247  (xin / xref  yin / yref)
    X and Y axis rotation: 359.99982  359.99982  (degrees  degrees)

Median sky level for mrgN20151214S0239:   2949.7
Transforming mrgN20151214S0240 to mrgN20151214S0240_trn

Transforming image mrgN20151214S0240[SCI] to image mrgN20151214S0240_trn
    Using transform mrgN20151214S0240_trn in database imcoadd.dat
Median sky level for mrgN20151214S0240:   2990.4
Transforming mrgN20151214S0241 to mrgN20151214S0241_trn

Transforming image mrgN20151214S0241[SCI] to image mrgN20151214S0241_trn
    Using transform mrgN20151214S0241_trn in database imcoadd.dat
Median sky level for mrgN20151214S0241:   2989.4
Transforming mrgN20151214S0242 to mrgN20151214S0242_trn

Transforming image mrgN20151214S0242[SCI] to image mrgN20151214S0242_trn
    Using transform mrgN20151214S0242_trn in database imcoadd.dat
Median sky level for mrgN20151214S0242:   2961.5
Median sky level for mrgN20151214S0239:   2949.7
Median sky level for mrgN20151214S0240:   2990.4
Transforming mrgN20151214S0239badpix +mask to mrgN20151214S0240badpix

Transforming image tmpmsk4365_1830 to image mrgN20151214S0240badpix
    Using transform mrgN20151214S0240_trn in database imcoadd.dat
Median sky level for mrgN20151214S0241:   2989.4
Transforming mrgN20151214S0239badpix +mask to mrgN20151214S0241badpix

Transforming image tmpmsk4365_1867 to image mrgN20151214S0241badpix
    Using transform mrgN20151214S0241_trn in database imcoadd.dat
Median sky level for mrgN20151214S0242:   2961.5
Transforming mrgN20151214S0239badpix +mask to mrgN20151214S0242badpix

Transforming image tmpmsk4365_1871 to image mrgN20151214S0242badpix
    Using transform mrgN20151214S0242_trn in database imcoadd.dat

Mean intensity for   mrgN20151214S0239_trn:  (N,absolute,relative)=   72  1198598  1.000
Mean intensity for   mrgN20151214S0240_trn:  (N,absolute,relative)=   72  1199217  1.001
Mean intensity for   mrgN20151214S0241_trn:  (N,absolute,relative)=   72  1201022  1.002
Mean intensity for   mrgN20151214S0242_trn:  (N,absolute,relative)=   72  1178268  0.983

Masking cosmic ray events in mrgN20151214S0239_trn
Applying growth radius 1
Masking cosmic ray events in mrgN20151214S0240_trn
Applying growth radius 1
Masking cosmic ray events in mrgN20151214S0241_trn
Applying growth radius 1
Masking cosmic ray events in mrgN20151214S0242_trn
Applying growth radius 1
IMCOADD exit status: good
-----------------------------------------------------------------------------

STEP 10: Inspect and save product


In [21]:
iraf.display(science['procsci']+'[sci,1]', 1, \
             zrange='no', zscale='no', z1=2000, z2=10000)

shutil.copy(science['procsci'], products_path)

# Final reduced image
os.listdir(products_path)


z1=2000. z2=10000.
Out[21]:
['ngc3310_g.fits', 'ngc3310_i.fits', 'ngc3310_r.fits']

In [ ]: